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We further study a predictive model for the masses and mixing matrix of three Majorana neu- 
trinos. At zeroth order the model yielded degenerate neutrinos and a generalized "tribimaximal" 
mixing matrix. At first order the mass splitting was incorporated and the tribimaximal mixing ma- 
trix emerged with very small corrections but with a zero value for the parameter S13. In the present 
paper a different, assumed weaker, perturbation is included which gives a non zero value for S13 and 
further corrections to other quantities. These corrections are worked out and their consequences 
discussed under the simplifying assumption that the conventional CP violation phase vanishes. It 
is shown that the existing measurements of the parameter S23 provide strong bounds on S13 in this 
model. 

PACS numbers: 14.60.Pq, 12.15.F, 13.10.-|-q 



I. INTRODUCTION 



At present, the particle physics community is planning, as a follow-up to the enormously important experiments of 
the last decade [H-0J extensive program with the goal of more accurately understanding the neutrino masses and 
mixings. There is really no accepted theory for an a priori prediction of these quantities. Hence it seems worthwhile 
to investigate in detail various theoretical models to develop plausible scenarios which might be tested. 

Here we look more closely at a particular model presented in Q and further studied in Q and in [l^l . That model 
assumed an initial permutation symmetry {S^) which is motivated by the fact that the 3x3 matrix which transforms 
the defining representation to irreducible form is, up to a single parameter rotation, the same as the "tribimaximal" 
matrix, which is in, at least, rough agreement with the present experimental situation. The tribimaximal form is 
taken to be: 
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The assumption was also made that, at zeroth order, the three neutrinos are degenerate. It may be seen from Table 
I of Q that this is plausible for a large range of possible fits to the data. However such an assumption at first seems 
inconsistent with permutation symmetry which suggests two of the three neutrinos to be degenerate (however not 
the two "solar" neutrinos) and different in mass from the third. The proposed solution to this problem called for 
the introduction of a Majorana type phase, which does not affect the usual neutrino oscillations but does affect the 
rate for neutrinoless double beta decay. The complications involved in obtaining a suitable Higgs scheme for both 
the neutrino mass matrix and the charged lepton mass matrix (which can be arranged to be proportional to the unit 
matrix) in this approach are discussed in some detail in 

Of course, many interesting different models for neutrinos based on permutation symmetry have been discussed 
for a long time [Tlj - [Tt} . In addition, many interesting models with similar approaches to the tribimaximal mixing 



matrix have been vigorously pursued [18 1 



In the model under present discussion, the zeroth order piece of the prediagonal Majorana neutrino mass matrix 
has the well known 5*3 invariant form: 
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Here a and /3 arCjin general, complex numbers while d is usually called the "democratic" matrix. As discussed in 
detail in [l^l and Q, we take 

a = -i\a\e-''^/^, (3) 

where the physical phase tp lies in the range: 

< -0 < 71". (4) 
For the assumed initial degeneracy, \a\ is related to f3, assumed real, by: 

3^ 



(5) 

The two zeroth order parameters are the degenerate neutrino masses, |a| and the phase tp which contributes to the 
neutrinoless double beta decay amplitude. 



2sm(V'/2)' 

The first order perturbation treated in [9|] and ^0|] is 




(6) 

where t and u are parameters. In general, t and u may be complex but they were assumed real for simplicity. This 
perturbation is well known as the "mu-tau" symmetry |3]|-(33j. The assumed zeroth order degeneracy (which actually 
may be relaxed, if desired) forces us to use degenerate perturbation theory. Then A turns out (eg, section II of [10]) 
to be the only possible choice which forces the desired tribimaximal form (as opposed to the generalized tribimaximal 
form) of the first order mixing matrix. 

Here we will choose for the second order perturbation, the matrix: 

/ <' u' 0\ 

A' = u' <' . (7) 

\ y 

For simplicity we again consider the parameters, t' and u' to be real. 

Note that this second order perturbation preserves the 5*2 subgroup which involves the 1-2 interchange. One might 
wonder about also including a perturbation. A" which preserves the 1-3 5*2 subgroup. However, that is not expected 
to give anything new since the combination of Eqs.([5]), ([B]) and ([7|) already has the same number of parameters as the 
most general symmetric matrix, M^. 

The combination of Eqs.(I2]), (O and (O was motivated by the group theory treatment of the strong interactions 
before QCD which led for example to the Gell-Mann Okubo mass formula |34]. In that case the initial term was 
flavor SU(3) invariant, the next term was invariant under the SU(2) isospin subgroup while the smallest last term 
was invariant under the SU(2) "U-spin" subgroup. In the present case the zeroth order term has the discrete group 
S3 invariance and two different S2 subgroups are left invariant by the two perturbations. 

II. PERTURBATION ANALYSIS 
In [To| we diagonalized the needed symmetric matrix: 

R^{al+ Pd + A)R = 

I t + u ^{t + u) \ 
al+ ^{t + u) 3/3+|(t + u) ■ (8) 
V t-uj 

The diagonalization of this matrix gave the first order neutrino mixing matrix, as 

X'l' = RRi, (9) 

where. 
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This results in a diagonalization with complex eigenvalues. To make these real positive we multiplied if*-^-* on the 
right by a suitable diagonal matrix of phases. 

To include the 2"-''-order perturbation, Eq. ([7]), we must diagonalize, 



H = R[R^{aI + I3d + A + A')RRi 



(11) 



where, after some computation and neglect of still higher order terms, we obtain: 
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We introduced the notation H'^ (Everything in Eq. pTj) except for A') and H' to indicate that, rather than making 
an explicit diagonalization we will regard, the result to first order as a "zeroth order Hamiltonian" , the given second 
order term, Eq. ([7]) as a "first order perturbation" and use ordinary quantum mechanics perturbation theory to proceed. 
In that approach one has of course the corrections to the energies as: 



(13) 



while the corrections to the eigenvectors are. 



(14) 



A more general perturbation approach, which gives the same results, is discussed in the Appendix. The lepton mixing 
matrix up to and including second order then reads: 



K = RRiR2P={iPi,iP2,^3)P, 



(15) 



where the ^jji are the columns of RR1R2 and furthermore P is the phase matrix needed for the neutrino masses to 
be real positive; explicitly. 
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and the phase matrix has the form. 
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wherein, 



TT 1 coUib/i) 



^ 9,3 90 90 

1 -ir coi(VV2) 1 
TT — —tan ' ' 



2 ^ 4_ 4(t+n) 4(t'+n') 

9/3 9/3 



TT 



3/3 3/3 

Note that we are free to subtract (t + cr + p)/3 from each of these three entries. Then the sum of the modified three 
entries will vanish in accordance with the requirement that there be only two independent Majorana phases. The real 
positive neutrino masses to second order are then: 



3 2, 5, , oV'i 

mi « -Pcsc-[l ~ — (t + u+-t -2u )sin -J, 

7712 ~ ^/3CSC^[1 + ^ + + ^' + 

W3 - l^P<^sc^[^ -—{t-u+ f')sin^^]. (19) 

Notice that the zeroth order masses have the characteristic strength, /? while the first order masses are suppressed 
by {t,u)/f3 and the second order masses are suppressed by {t',u')/f3. 

Also notice that the absolute values of the neutrino masses depend on the Majorana phase, -0. However, the lepton 
number conserving neutrino oscillations can not depend on a Majorana phase [Sa]- As a check of this we see that the 
phase ip cancels out when one considers the mass differences ^ 

g 

A = ml- ml Pi 3/3(t + m) + -fit', 

7 

B = 771^ - 777^ « I3{~bt + u)~ l3{-t' + 2u'), 

C = ml-ml^ 2P{-t + 2u) + /3(t' - 2u'). (20) 
Of course, A, B and C are not independent. There are two, presently unresolved, experimental possibilites: 



Typel : 
Type2 : 



7773 > 7772 > 777 1, 
7772 > ^771 > "7,3. 



(21) 



The corresponding relations are: 



Typel: \C\ = \B\ + A, 

Type2 : \C\ = \B\ - A. 



(22) 



These relations were obtained by using the known positive sign of A and that only the two possibilities 7773 > 7773 > 
mf and 7772 > 777^ > 777I are allowed. In the literature some works specify A and \B\ while others specify A and \C\. 
The following best fit values for the perturbation parameters /3t and /3u were given in the first order treatment [lOl]: 
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(3t « -4.13 X IQ-^eV'^, 

Pu K. 4.39 X IQ-^eV"^, Typel (23) 



/3t w 4.21 X lO^^eV^, 

I3u w -3.94 X IQ-^eV'^ Type2. (24) 

III. ELEMENTS OF THE MIXING MATRIX 

We employ the following parameterization [36] of the leptonic mixing matrix, K: 

C12C13 S12C13 

K = \ -S12C23 - Ci2Si3S23e*^ C12C23 - Si2Sl3S23e''' C13S23 | -P, (25) 
S12S23 - Cl2Sl3C23e*^ -C12S23 

where C12 is short for cos9i2 for example. P is the diagonal matrix of Majorana type phases given in Eqs. (jl7p and 
([TS]) for the present model. For simplicity we are presently neglecting the conventional CP violation and thus setting 
7 = 0. To specify S12, S13 and S23, it is clearly sufficient to compare the (1-2), (1-3) and (2-3) matrix elements of K 
in Eq. ips)) with those calculated in Eg. p^ . This yields: 

1 2 t + u 1 t' + u' 

S12C13 
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For an initial orientation we see that at zeroth order, S13 vanishes and also K has the tribimaximal form. When 
the first order perturbation characterized by t and u is added, neither S13 nor S23 change. However S12 is somewhat 
modified as discussed previously in section IV of [loj . When the second order perturbation characterized by t' and u' 
is added, S13 finally becomes non-zero while both S12 and S23 suffer further corrections. 

But something unusual is happening; there are terms for S13 and S23 which behave like t' /t and are manifestly of 
first order in strength. These arise from the energy difference denominator in Eq. (|14p . Since we had to use degenerate 
perturbation theory at first order this denominator is proportional to the first order "energy" corrections rather than 
the zeroth order energies. Keeping terms of actual first order in strength we find the interesting relation; 

si3 « -2,5s23, (27) 
where SS23 denotes the deviation of S23 from its tribimaximal value. Also the good approximation C13 = 1 was made. 



IV. NUMERICAL ESTIMATES 

Already, Fogli et. al. [13] and Schwetz et. al. [l^ have pointed out that detailed analysis of existing neutrino 
oscillation experiments gives some hint for non zero S13. Thus it seems interesting to see what predictions emerge 
from Ea.(P7|. 

Expanding S23 around its "tribimaximal value" as S23 — [s23]tbm + ^S23, one gets: 

(S23)' + V26s23- (28) 

Comparing with the results of a global analysis of the oscillation data given in Table Al of [s^ one then identifies, 
for respectively la, 2a and 3a errors: 

)5s23|= 0.05, 0.08, 0.11. (29) 
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Note that the three cases are associated with the experimental data relating to the 2-3 type neutrino oscillations. 
Using Eq. (P7|) then leads to the corresponding predictions, 

|si3| < 0.025, 0.040, 0.055. (30) 

It is amusing to note that these values range from about 1/4 to 1/2 of the "best fit" value jsiaj = 0.11, which is also 
presented in the first column of Table Al in [s^. Of course, our estimates provide a test of the present theoretical 
model for neutrino parameters and have no connection with experimental data on |si3|. 

As discussed above, the theoretical estimate for |si3|, is of characteristic first order strength, appearing as a ratio 
of a second order quantity divided by a first order quantity. Using Eg. ([^5]) for sia and neglecting the term of second 
order strength we can get an estimate of the relative second to first order effects: 

lV^I~2V2|si3|« 0.071, 0.11, 0.16, (31) 
t — 2u 

wherein Eg. ([50]) was used. Evidently the second order effects seem to be suppressed by about 1/10 compared to 
the first order effects. On the other hand, as seen in Eg. ipOl) . the quantities t' and u' enter in the true second order 
corrections for the neutrino mass differences. Thus those corrections are likely to be small- on the order of ten percent 
of the first order mass splittings. 



V. SUMMARY AND DISCUSSION 



In this work, we designated the zeroth order parameter as /3, the first order parameters as t and u and the second 
order parameters as t' and u' . The first order corrections to the neutrino masses were suppressed by (t, u)/ j3 compared 
to zeroth order. For the mixing angles, the first order corrections had a previously obtained piece proportional to 
{t,u)/f3 as well as a new piece proportional to {t' ,u')/ {t,u). The latter term arose because we are using degenerate 
perturbation theory and is clearly important for sia to be non-zero and correlated to corrections of S23- 

Here, we have numerically neglected, for both masses and mixing angles terms proportional to {t',u')/l3. In [To| 
we considered {t,u)/l3 to be about 1/5. Here we found a characteristic strength of S13 to correspond to {t' ,u') /{t,u) 
about 1/10. Both of these magnitudes are roughly similar. 

Note that Egs. (PD|) for the neutrino mass differences and Egs. (P51) for the mixing angles do contain pieces of actual 
second order strength. These should be interesting to study in the future when more precise data becomes available. 

The first order corrected formula for the neutrinoless double beta decay factor is given in Eq.(51) of [3. This was 
derived from Eq.(49) in which (sia)^ was set to zero. Now S13 is not zero but its square contributes at a higher order. 
Furthermore, it is easy to see, using Egs. p^ . that the first two terms in Eq.(49) do not have any contributions of 
first order strength like {t' ,u')/ {t,u). Hence that formula for rriee still holds to first order. 
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Appendix A: Alternative perturbation method 



We present here an alternative approach which leads to results in perturbation theory order by order. This can be 
applied to the case at hand or more generally when the mass matrix is invariant at zeroth order under a finite group 
Go and then we add perturbations of decreasing importance in the small parameter x such that for example the n*'* 
perturbation is of order x" and is invariant under a smaller group Gn- The mass matrix can then be written as an 
expansion in x, 

M{x) = Mo + xMi + x^M2 + ... (Al) 
where Mq is invariant under Go, Mi under Gi and so on. 
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The eigenvalues (diagonal) and eigenvector matrices can also be expanded as, 

Md{x) = Mdo + xMdi + x'^Md2 + .: 

R{x) = Ro + xRi+x^R2 + ... (A2) 

where, 

R^{x)M{x)R{x) ^ Md{x) (A3) 

is the eigenvalue equation. 

If we differentiate Eq (jA3[) once we obtain: 



R^'MR + R^M'R + R^MR' = M'^ (A4) 

which can be written as: 

[Md, R^R'] + R^M'R = M'd (A5) 

Here we used the orthonormality condition for the eigenvector matrix: 

R^'R + R^R' = (A6) 

Note that the matrix R^'R which appears in what follows is antisymmetric (in each order of perturbation theory) 
and in consequence all of its derivatives will be antisymmetric. 
The second derivative and third derivative equations will read: 

[Md, R^R'] + [Md, {R^R')'] + [R^M'R, R^R'] + R^ M" R = M'J, 
[M^, R^R'] + 2[M'd, {R^R')'] + [M<i, {R^R')"] + 

[[R^M'R, R^R'],R^R'] + 2[R^M"R, R^ R'] + [R^M'R, [R^ R')'] + R^M"'R = M'J' (A7) 

All commutators of diagonal matrices give zero on diagonal and in consequence the mass eigenvalues are obtained 
from the rest of the terms. 

It is clear that by setting x = one can associate the first derivative with the first order perturbation theory, 
second with second order and so on. The mass eigenvalues and the matrix R^ R' can be extracted in each order from 
equations like Eq (IX5|) and Eq (|X7)) . 

Then one should use the orthonormality condition to obtain the eigenvector matrix according to: 

R^{x)R'{x) = RIRi + x{rIRi + 2i?Ji?2) + .... (A8) 

Using this method and Gq = 5*3, Gi = 82^ and G2 — S12 one retrieves the eigenvalues and eigenvectors in each 
order of perturbation theory. The results agree with those presented in the main text. 
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